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ABSTRACT 

We report results from the first deep millimeter continuum survey targeting Brown Dwarfs (BDs). The survey 
led to the first detection of cold dust in the disks around two young BDs (CFHT-BD-Tau 4 and IC 348 613), with 
deep JCMT and IRAM observations reaching flux levels of a few mJy. The dust masses are estimated to be a few 
Earth masses assuming the same dust opacities as usually applied to TTauri stars. 

Subject headings: accretion, accretion disks — circumstellar matter — stars: formation - stars: low-mass, brown 
dwarfs - stars: individual (CFHT-BD-Tau 4, IC 348 6 1 3) 



1. INTRODUCTION 

Several near-infrared surveys (Muench et al. 2001; Oliveira 
et al. 2002; Liu et al. 2003; Jayawardhana et al. 2003) and mid- 
infrared measurements (Comeron et al. 2000; Persi et al. 2000; 
Testi et al. 2002; Natta et al. 2002; Apai et al. 2002) indicate 
the presence of disks around Brown Dwarfs (BDs). In contrast 
to infrared emission, submillimeter and millimeter emission is 
certainly always optically thin and is an excellent measure for 
the total dust mass. 

An early attempt to detect millimetre continuum emission 
also from very low-mass stars and suspected BDs was done 
by Andre & Montmerle (1994). Probably the most sensitive 
observation of this kind was carried out by Carpenter (2002) 
using the OVRO interferometer, however at a relatively long 
wavelength. We carried out the first successful search for dust 
continuum emission associated with confirmed BDs, using the 
bolometer arrays SCUBA at the JCMT and MAMBO at the 
IRAM 30-m telescope. 

The survey led to the detection of circumstellar dust around 
the two young BDs CFHT-BD-Tau 4 and IC 348613, which 
have ages below 10 Myrs. In the case of field BDs we obtained 
upper mass limits of a few Moon masses of dust. For BDs in the 
Pleiades the mass limits are less strict and range 4 and 7 Earth 
masses. We should note that the data presented for CFHT-BD- 
Tau 4 together with other ground-based and ISO data allowed 
the first detailed discussion of a complete spectral energy distri- 
bution of a BD, ranging from optical to millimeter wavelengths 
(Pascucci et al. 2003). 

The detection of these amounts of circumstellar material 
around two young BDs makes the formation of planets or even 
planetary systems around BDs a possibility. Therefore, search 
strategies for planets should include BDs. In the case of imag- 
ing surveys, they might even be the best targets. 

In addition, the detection of significant amount of circum- 



stellar dust carries important information about the formation 
processes of BDs. These detections, together with the discov- 
ery of quite a number of BD binaries (Bouy et al. 2003; Bur- 
gasser et al. 2003), do not support fragmentation of circumstel- 
lar disks as the general process for BD formation. Other for- 
mation scenarios for BDs (Bate et al. 2003; Reipurth & Clarke 
2001; Watkins et al. 1998) include ejection from multiple sys- 
tems and erosion of star-forming cloudlets by stellar winds and 
UV radiation from massive stars. Disks will certainly have dif- 
ferent structures depending on the formation mechanism. How- 
ever, statistics is still poor and information about the disk struc- 
ture from interferometric observations are needed before one 
can put more definite observational constraints on BD forma- 
tion scenarios. 

2. OBSERVATIONS 

To search for circumstellar material around BDs, we selected 
relatively young BDs with an age of a few Myrs because these 
objects should have the highest probability to be associated with 
disk material. The nine selected objects are located in Taurus, 
the a Orionis cluster, IC 348, and the Upper Scorpius OB asso- 
ciation (Martin et al. 2001; Bejar et al. 2001; Najita et al. 2000; 
Luhman 1999; Ardila et al. 2000). For the first three regions, 
additional selection criteria were the previous detection of Ha 
emission (Martin et al. 2001), the presence of X-ray emission 
(Mokler & Stelzer 2002; Preibisch & Zinnecker 2001), and the 
requirement that the objects should be as isolated as possible 
in order to avoid confusion during the observations with single- 
dish telescopes. 

The three BDs in the Taurus star-forming region are among 
the youngest BDs of our target list. They have ages of about 
1 Myr (Martin et al. 2001). The object CFHT-BD-Tau 4 shows 
the highest Ha emission among the Taurus BDs (Martin et al. 
2001), emits X-ray radiation (Mokler & Stelzer 2002) and 
shows mid-infrared excess emission (Pascucci et al. 2003). 
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The second group, we selected, are objects in the Pleiades 
which have a distance of 1 16 pc and an age of 0. 12 Gyr. The ob- 
jects are considerably older than the first group of BDs. In addi- 
tion, we searched for dust emission from very nearby field BDs. 
On the average, these objects are older than the first group, too, 
but the actual age determination contains large uncertainties. 

The targets with their ages and distances are compiled in Ta- 
ble 1. 

The observations were carried out at the James Clerk 
Maxwell Telescope (JCMT) on Mauna Kea, Hawai'i, and at 
the IRAM 30 m telescope, Pico Veleta, Spain. At the JCMT we 
used the 37-channel array camera SCUBA (Holland et al. 1999) 
at a wavelength of 850 /im. At the IRAM 30-m telescope, we 
performed the observations at 1.3 mm with the 37-channel ar- 
ray camera MAMBO (Kreysa et al. 1998). These wavelengths 
are a good tradeoff between the expected flux and the opacity 
of the atmosphere. 

The observations were performed between March 2002 and 
January 2003. The chosen observing mode is photometry 
with the central bolometer of SCUBA and MAMBO, respec- 
tively. The background (telluric and astronomical) subtraction 
is achieved by chopping with the secondary mirror and nodding 
the telescope. Its level has been estimated from the inner ring 
consisting of six bolometer pixels. 

The measured flux densities and derived upper limits are 
listed in Table 2. For many sources, we can only derive upper 
limits since the signal-to-noise-ratio (SNR) is less than three. 
The upper limits are given as 3a values, where a is the stan- 
dard deviation of the statistical error. 

3. RESULTS 

We detected millimeter continuum emission associated with 
two BDs. Both, BD Tau4 and IC 348 613, were detected at two 
wavelengths, i.e. with the JCMT at 850 fim and with the IRAM 
30m telescope at 1.3 mm. These detections can be used to esti- 
mate the amount of circumstellar matter around these BDs for 
the first time. 

Carpenter (2002) however using the OVRO interferometer to 
map IC 348 did not detect IC 348613. The reason is that he 
mapped at a wavelength of 3 mm. If we extrapolate the mea- 
sured flux densities to a wavelength of 3 mm, the flux density 
falls just below the detection limit of the OVRO map. 

3.1. Dust Masses 

The measured flux density F v at submillimeter/millimeter 
wavelengths is certainly optically thin thermal emission by cold 
dust heated by the BD. The dust mass M can be derived using 
the formula 



B v {T,\)n v {\y 

where F v stands for the flux density, B u for the Planck function, 
T for the dust temperature and k v for the mass absorption co- 
efficient. The quantity D is the distance to the object. Likewise 
upper limits for the flux density are translated into upper limits 
for the amount of circumstellar matter. 

Among the above quantities, the distances are well estab- 
lished from trigonometric parallaxes in the case of most of the 
field BDs. For the young BDs we used the known distances to 
the clusters/star-forming regions. However, the dust properties 
k v and the dust temperature T need a more thorough discus- 
sion. 



The targeted BDs have a large spread in age and we would 
expect different evolutionary stages of the disk material if 
present at all. Therefore, we will apply two different sets of 
dust parameters (see, e.g., Henning et al. 1995). The "young" 
dust parameters will be applied to BDs with ages up to 10 Myr. 
The "debris" dust parameters will be used for the other objects, 
which are all older than 100 Myr. We assume a constant dust 
temperature as an approximation to the temperature distribution 
of the circumstellar dust. The next paragraphs discuss the two 
sets of dust parameters. 

"Young" dust: For dust around young BDs, we choose, 
for the sake of comparison, a mass absorption coefficient of 
k v = 2cm 2 g _1 at 1.3 mm. The same value of k v and a gas-to- 
dust ratio of 1 00 was applied by Beckwith et al. ( 1 990) to derive 
disk masses for T Tauri disks. For the measurements at 850 /im, 
we assume a wavelength dependence of k v oc \~P with (3 = 1, 
also in accordance with Beckwith et al. (1990). This leads to 
k v = 3 cm 2 g _1 at 850 /zm. 

The plausible range of dust temperatures is relatively small. 
We assume an average temperature of 10 to 20 K for the dust. 
This is the range of the mass-averaged dust temperature in the 
models for the disk around CFHT-BD-Tau 4 discussed by Pas- 
cucci et al. (2003). 

"Debris" dust: The BDs in the Pleiades and the field objects 
are presumably older than 100 Myr. At this age one can no 
longer assume a TTauri-like disk with its dust properties. If 
there is a disk left it will resemble debris disks like the (3 Pictoris 
disk. To allow comparisons to submillimeter observations of 
debris disks, we adopted the mass absorption coefficient used 
by Dent et al. (2000), i.e. k v = 0.4. . . 1.7cm 2 g _1 at 850 /im. 
The gas-to-dust mass ratio in debris disks is controversially 
discussed (Thi et al. 2001; Lecavelier des Etangs et al. 2001). 
Therefore, we will only give the dust masses in these cases. 

Since debris disks are optically thin, the circumstellar dust 
grains around the old BDs have at least the temperature of in- 
terstellar grains, i.e. 20 to 30 K depending on the composition 
and the actual interstellar radiation field at the position of the 
BDs. The heating of the dust by the old BD is negligible, thus, 
we assume a temperature between 20 to 30 K. 

Applying the above discussed dust properties to the millime- 
ter continuum measurements, we obtain the dust masses com- 
piled in Table 2. The ranges result from the uncertainties in the 
mass absorption coefficient and the temperature. Furthermore, 
the millimeter emission has been measured at two wavelengths 
for the detections. The mass estimates derived from the two 
measurements differ slightly. For the detections, the masses in 
Table 2 correspond to the minimum and maximum mass esti- 
mates. In the case of non-detections, the mass limits are derived 
from the 3cr flux limits using the dust parameters yielding the 
highest masses. 

The mass of circumstellar dust for CFHT-BD-Tau 4 is be- 
tween 1.4 and 7.6 M E (Earth mass) and for IC 348 613 the mass 
ranges between 5.4 and 18 M E . The dust is certainly distributed 
in the form of disks as a detailed analyses of CFHT-BD-Tau 
4 by Pascucci et al. (2003) show. The disk masses for the two 
BDs are 0.4 . . . 2.4Mj (Jupiter mass) and 1 .7 . . . 5.7Mj, if we ex- 
trapolate the dust masses to disk masses assuming a gas-to-dust 
ratio of 100. The upper limits on the circumstellar dust mass 
around the other young BDs are around a few Earth masses. 

The dust mass constraints for the nearby field BDs are even 
more stringent, reaching upper limits of some Moon masses. 
Thus, even around the low-luminosity BDs circumstellar disks 
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cannot be long-lived. Debris disks around the old BDs, which 
may be associated with planetary systems, cannot be ruled out, 
since we were not sensitive down to fractions of Moon masses 
of dust which is the order of the dust mass around the low-mass 
main-sequence star e Eridani (Dent et al. 2000). 

3.2. Disks aroundyoung BDs 

One important question about BDs and their disks is the mass 
ratio of the central objects and the circumstellar material. The 
BDs' masses themselves (CFHT-BD-Tau 4: M* = 30...75M, 
\ Martin et al. 2001; IC 348 613: M* = 20. . .40Mj, Preibisch 
& Zinnecker 2001) are uncertain as well as the disk masses. 
Therefore we do not attempt to estimate their ratio here. Still, 
we note that the observations are consistent with a disk/BD 
mass ratio of a few percent as it is found for most of the T Tauri 
stars. 

The detection of some Jupiter masses of matter around BDs 
is especially interesting for extrasolar planet searches. Guen- 
ther & Wuchterl (2003) mention that planets may form around 
a BD if a disk of enough mass is present. Their radial velocity 
survey targets very low-mass stars and BDs because these tar- 
gets are relatively inactive. Furthermore, young planets around 
BDs would be the easiest to detect by direct imaging because 
the contrast between the BD and its planet would be much more 
favorable than for a planet around a higher-mass star. 

4. SUMMARY 

The observing campaign targeting BDs of several popula- 
tions resulted in detections of millimeter emission associated 
with two young BDs. The mass estimates for the young BD 



disks are 0.4...2.4M, and 1.7...5.7M, for CFHT-BD-Tau 4 
and IC 348 613, respectively. For the other targets, stringent up- 
per limits on the amount of circumstellar matter were derived 
from the measured upper limits on the millimeter continuum 
flux densities. To estimate the dust masses, two sets of dust 
properties had to be applied: "Young" dust properties to BDs 
younger than lOMyr and "debris" dust properties to BDs older 
thanlOOMyr. 

The detection of a few Jupiter masses of circumstellar matter 
around young BDs is an important result. To ensure this mass 
estimate, the dust properties have to be constrained further. 
However, a refinement of the dust properties will hardly change 
the fact that there are substantial amounts of circumstellar mate- 
rial around the two BDs CFHT-BD-Tau 4 and IC 348 613. Thus, 
the detections make BDs to places of possible planet forma- 
tion. This fact opens a new set of targets for extrasolar planet 
searches, especially for direct imaging because of the low con- 
trast between the central object and an prospective planet. 
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acknowledges support through the DFG grant He 1935/15-1. 



REFERENCES 



Andre, P. & Montmerle, T. 1994, ApJ, 420, 837 

Apai, D., Pascucci, I., Henning, T., Sterzik, M. F., Klein, R., Semenov, D., 

Gunther, E., & Stecklum, B. 2002, ApJ, 573, LI 15 
Ardila, D., Martin, E., & Basri, G. 2000, AJ, 120, 479 
Bejar, V. J. S. et al. 2001, ApJ, 556, 830 

Bate, M. R., Bonnell, I. A., & Bromm, V. 2003, MNRAS, 339, 577 
Beckwith, S. V. W., Sargent, A. I., Chini, R. S., & Guesten, R. 1990, AJ, 99, 
924 

Bouy, H., Brandner, W., Martin, E. L., Delfosse, X., Allard, E, & Basri, G. 

2003, AJ, submitted (astro-ph/0305484) 
Burgasser, A. J., Kirkpatrick, J. D., Reid, I. N., Brown, M. E., Miskey, C. L., & 

Gizis, J. E. 2003, ApJ, 586, 512 
Carpenter, J. M. 2002, AJ, 124, 1593 

Comeron, E, Neuhauser, R., & Kaas, A. A. 2000, A&A, 359, 269 
Dahn, C. C. et al. 2002, AJ, 124, 1 170 

de Zeeuw, P. T., Hoogerwerf, R., de Bruijne, J. H. J., Brown, A. G. A., & 

Blaauw, A. 1999, AJ, 1 17, 354 
Dent, W. R. F, Walker, H. J., Holland, W. S., & Greaves, J. S. 2000, MNRAS, 

314, 702 
Festin, L. 1998a, A&A, 333, 497 
— . 1998b, MNRAS, 298, L34 

Freed, M., Close, L. M., & Siegler, N. 2003, ApJ, 584, 453 

Gliese, W. & JahreiB, H. 1991, Preliminary Version of the Third Catalogue of 

Nearby Stars, Tech. rep. 
Guenther, E. W. & Wuchterl, G. 2003, A&A, 401, 677 

Henning, T, Michel, B., & Stognienko, R. 1995, Planet. Space Sci., 43, 1333 
Holland, W. S., Robson, E. I., Gear, W. K., Cunningham, C. R., Lightfoot, J. F, 
Jenness, T., Ivison, R. J., Stevens, J. A., Ade, P. A. R., Griffin, M. J., Duncan, 
W. D., Murphy, J. A., & Naylor, D. A. 1999, MNRAS, 303, 659 
Jayawardhana, R., Ardila, D. R., Stelzer, B., & Haisch Jr, K. E. 2003, aJ, in 
press 

Kirkpatrick, J. D. et al. 2000, AJ, 120, 447 

Kreysa, E. et al. 1998, in Proc. SPIE, Vol. 3357, 319 

Lecavelier des Etangs, A. et al. 2001, Nature, 412, 706 

Leggett, S. K., Allard, F, & Hauschildt, P. H. 1998, ApJ, 509, 836 



Leggett, S. K., Toomey, D. W., Geballe, T. R., & Brown, R. H. 1999, ApJ, 517, 
L139 

Liu, M. C, Najita, J., & Tokunaga, A. T. 2003, ApJ, 585, 372 
Luhman, K. L. 1999, ApJ, 525, 466 

Martin, E. L., Brandner, W., & Basri, G. 1999, Science, 283, 1718 

Martin, E. L., Dougados, C, Magnier, E., Menard, F, Magazzu, A., Cuillandre, 
J.-C, & Delfosse, X. 2001, ApJ, 561, L195 

Mermilliod, J.-C, Turon, C, Robichon, N., Arenou, E, & Lebreton, Y. 1997, 
Proceedings of the ESA Symposium 'Hipparcos - Venice '97', 13-16 May, 
Venice, Italy, ESA SP-402 (July 1997), p. 643-650, 402, 643 

Mokler, F. & Stelzer, B. 2002, A&A, 391, 1025 

Muench, A. A., Alves, J., Lada, C. J., & Lada, E. A. 2001, ApJ, 558, L51 

Najita, J. R., Tiede, G. P., & Carr, J. S. 2000, ApJ, 541, 977 

Nakajima, T, Oppenheimer, B. R., Kulkarni, S. R., Golimowski, D. A., 

Matthews, K, & Durrance, S. T. 1995, Nature, 378, 463 
Natta, A., Testi, L., Comeron, E, Oliva, E., D'Antona, F, Baffa, C, Comoretto, 

G., & Gennari, S. 2002, A&A, 393, 597 
Oliveira, J. M., Jeffries, R. D., Kenyon, M. J., Thompson, S. A., & Naylor, T. 

2002, A&A, 382, L22 
Oppenheimer, B. R., Kulkarni, S. R., Matthews, K, & Nakajima, T. 1995, 

Science, 270, 1478 

Pascucci, I., Apai, D., Henning, T, & Dullemond, C. P. 2003, ApJ, 590, LI 11 
Persi, P. et al. 2000, A&A, 357, 219 

Preibisch, T, Guenther, E., Zinnecker, H, Sterzik, M., Frink, S., & Roeser, S. 

1998, A&A, 333, 619 
Preibisch, T. & Zinnecker, H. 2001, AJ, 122, 866 
Reipurth, B. & Clarke, C. 2001, AJ, 122, 432 
Ruiz, M. T., Leggett, S. K., & Allard, F. 1997, ApJ, 491, L107 
Scholz, R.-D. et al. 1999, A&AS, 137, 305 

Stauffer, J. R., Schultz, G., & Kirkpatrick, J. D. 1998, ApJ, 499, L199 

Testi, L., Natta, A., Oliva, E., D'Antona, F, Comeron, E, Baffa, C, Comoretto, 

G., & Gennari, S. 2002, ApJ, 571, L155 
Thi, W. F. et al. 2001, Nature, 409, 60 

Watkins, S. J., Bhattal, A. S., Boffin, H. M. J., Francis, N., & Whitworth, A. P. 
1998, MNRAS, 300, 1214 



1 Martin et al. (2001) derives a mass of 0.03 M© for CFHT-BD-Tau 2 and 3 with spectral types M8 and M9, respectively. The CFHT-BD-Tau 4 has the spectral type 
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Table 1 



Target list 



No. 


Target 




Co-ordinates 


Dist. 


Age 


References 






RA(2000) 


DEC(2000) 


(pc) 


(Myr) 






BDs in Taurus 








140 


1 


14 


1... 


CFHT-BD-Tau 1 


04*34" 


15.2 s 


+22°50'31" 








2... 


CFHT-BD-Tau 2 


04*36" 


10.4 s 


+22°59'56" 








3... 


CFHT-BD-Tau 4 
BD in <t Ori cluster 


04*39" 


47.3 s 


+26°01'39" 


370 


1 


17 

2, 15 


4... 


S Ori 03 
BDinIC348 


05*39" 


20.8 s 


-02° 30' 35" 


260 


0.5-10 


10, 21 


5... 


IC 348 613 

BDs in U. Scorpio. 


03*44" 


26.9 s 


+32°09'24.8" 


145 


5 


1, 18, 22 


6... 


USco 100 


16*02" 


04.13 s 


-20°50'41.5" 








7... 


USco 104 


15*57" 


12.66 s 


-23°43'45.3" 








8... 


USco 112 


16*00™ 26.57 s 


-20°56'32.0" 








9... 


USco 128 


15*59' 


11.20 s 


-23°37'59.0" 










BDs in the Pleiades 








116 


120 


4, 13, 20 


10... 


NPL36 


03*48" 


19.1 s 


+24° 25' 15" 








11... 


NPL37 


03*47" 


12.1 s 


+24°28'31" 








12... 


NPL38 


03*47' 


50.4 s 


+23°54'49" 








13... 


NPL 30 (Teide 1) 
Field BDs 


03*47' 


17.9 s 


+24°22'32" 




>100 a 




14... 


Gl 229B 


06*10" 


35.1 s 


-21°51'18" 


5.8 


500 


9, 12, 16 


15... 


2MASSI J0746425+200032 b 


07*46" 


42.5 s 


+20°00'32" 


12.3 




3,7 


16... 


2MASS1 J0825 196+21 1552 


08*25" 


19.6 s 


+21°15'52" 


12.5 




3,7 


17... 


LHS 2397a"' 


11*21' 


49.3 s 


-13°13'09" 


14.2 


2-12 Gyr 


5, 6, 7, 8 


18... 


Kelu 1 


13*05" 


40.2 s 


-25°41'06" 


19.2 


0.3-1 Gyr 


3, 7, 11, 19 


19... 


TVLM 868-110639" 


15*10' 


17.2 s 


-02°41'07" 


25.4 




3, 7,8 





"Lower limit for the objects where no individual age estimate is available. Comparison with evolutionary tracks 
suggests ages of the order of 100 Myr and higher, see (3). 

b Binary BD 

L BD Companion (15) 

"Objects at the substellar boundary 
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Table 2 



Observational results 





PI liv Cfl) 1 ^ m m 
nilAV"' 1 .J 111111 


Dust mass 




yiiuy ) 


(M E ) 


BDs in Taurus 








1... 


<2.73 


< 


4.8 


2... <5.81 




< 


4.1 


3... 10.8 ±1.8 


2.1±0.6 


1.4. 


..7.6 


BD in a Ori cluster 








4... 


<2.97 


<32 


BD in IC 348 








5... 7.6 ±2.4 


2.8±0.8 


5.4. 


.18 


BDs in Upper Scorpio 








6. . . <4.09 




< 


3.1 


7... <5.27 




< 


3.8 


8... <6.11 




< 


4.4 


9... <8.35 




< 


6.3 


BDs in the Pleiades 








10... <6.04 




< 


6.7 


11... <4.18 




< 


4.6 


12. . . <4.04 




< 


4.4 


13... <4.39 




< 


4.9 


Field BDs 








14... <9.41 




< 


0.03 


15... <3.76 




< 


0.05 


16... <5.16 




< 


0.07 


17... <4.79 




< 


0.08 


18... <2.15 




< 


0.06 


19... <6.73 




< 


0.36 



